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ABSTRACT 

We study a model of fermions interacting with a gauge field and calculate 
gauge-invariant two-particle Green's functions or response functions. The leading 
singular contributions from the self-energy correction are found to be cancelled 
by those from the vertex correction for small q and O. As a result, the remaining 
contributions are not singular enough to change the leading order results of the 
random phase approximation. It is also shown that the gauge field propagator is 
not renormalized up to two-loop order. We examine the resulting gauge-invariant 
two-particle Green's functions for small q and but for all ratios of Q/vFq 
and we conclude that they can be described by Fermi liquid forms without any 
diverging effective mass. 
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I. INTRODUCTION 

The problem of two dimensional fermions coupled to a gauge field has been a recent 
subject of intensive research. This problem appears as a low energy effective model of two 
different strongly correlated electronic systems, i.e., electrons in the fractional quantum 
Hall (FQH) regime and the high-temperature superconductors (HTSC), both of which 
have been considered as one of the most important problems in modern condensed matter 
physics. 

As the first example, this problem arises in a theory of the half-filled Landau level 
(HFLL) [1-3] in connection with the composite fermion theory of the FQH effect [4]. A 
composite fermion is generated by attaching even number of flux quanta to an electron 
[4]. The transformation from the electron to the composite fermion can be realized by 
introducing an appropriate Chern-Simons gauge field [1,5]. Especially, at the filling fraction 
v = 1/2, composite fermions see effectively zero magnetic field at the mean field level [1-4] 
because of the cancellation between the average of the Chern-Simons gauge field (from the 
attached magnetic flux quanta) and the external magnetic field. Thus, at the mean field 
level the system can be described as a Fermi liquid of composite fermions. The fluctuation 
of the gauge field beyond the mean field level has been studied within the random phase 
approximation (RPA) [1,3], which explains qualitative features of the recent experiments 
[6-11]. 

The other source comes from the recent gauge theory of the normal state of high 
temperature superconductors [12-15]. The gauge field arises as a fluctuation of the spin 
chirality [12] above the uniform resonating-valence-bond mean field state of the t — J 
model which is supposed to be an effective model of HTSC. The origin of the gauge field 
fluctuation can be traced back to the constraint that the doubly-occupied sites are not 
allowed because of the strong on-site Coulomb repulsion [12,13]. It has been suggested 
that the gauge field fluctuations play important roles in explaining anomalous transport 
properties of the normal state of HTSC [12,15,16]. 



Besides these real examples, the problem of fermions interacting with a gauge field has 
been studied as a potential example of non- Fermi liquids [17-28]. In contrast to the usual 
long-range interactions such as the Coulomb interaction, the transverse part of the gauge 
field cannot be screened because the gauge invariance requires the gauge field to be massless 
in the absence of symmetry breaking [17-19]. Thus, one can expect that the long-range 
interaction due to the transverse part of the gauge field gives rise to non-Fermi-liquid-like 
behaviors. In fact, some singular behaviors appear in the lowest order self-energy correction 
of fermions by the gauge field fluctuation [14,17-20]. The singular self-energy correction 
makes perturbative calculation unreliable at low energies. This motivated several non- 
perturbative calculations of one-particle Green's function of fermions which show highly 
non-Fermi-liquid-like behaviors [21-24,26]. It turns out that, even in the lowest order, the 
singular self energy correction makes the effective mass of the fermion divergent so that 
the usual single particle picture breaks down [1]. 

However, recent experiments on the electrons in the half-filled Landau level showed 
essentially Fermi-liquid-like behaviors [6-11] and also measured finite effective mass of 
composite fermions [10]. Therefore, we are in a situation that experiments apparently 
contradict to the insight we got from the one-particle Green's function of the fermions. 
However, the one-particle Green's function for the fermions is not gauge invariant. The 
singular self-energy correction in the one-particle Green's function (which leads to divergent 
effective mass [1]) may be an artifact of the gauge choice rather than a property of physical 
quasi-particles. Since it is not a gauge-invariant quantity, the one-particle Green's function 
for the fermions cannot be directly measured in experiments. It is possible that some 
singularities in the gauge-dependent one-particle Green's function simply do not appear 
in gauge-invariant correlation functions. One purpose of this paper is to examine some 
gauge-invariant response functions in order to determine whether the singular behaviors in 
the one-particle Green's function appear in gauge-invariant correlation functions or not. 

The importance of the gauge-invariance in calculating correlation functions can be 
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also seen in the following example. The leading order corrections (two-loop order) to the 
transverse polarization function (or current- current correlation function) are given by the 
diagrams in Fig. 3. Note that the sum of contributions from Fig. 3 (a)-(d) is not gauge- 
invariant because they contain only self-energy corrections but do not contain the vertex 
correction. For concreteness, let us consider the case of r\ = 2 in the model given by Eq.(8), 
which corresponds to the case of HTSC and the short-range interaction between fermions 
in HFLL. We also consider O <C vpq and q <^kp limits. In this case, it can be shown that 
the correction to the transverse polarization function due to the self-energy corrections 
(given by Fig. 3 (a)-(d)) has the following form: 



„ mn . l(q ,n)^A(2», (1) 

27T7 VFq kpq 



while the contribution from the free fermions is given by 

lmn? l(q ,0) = -^ SL (2) 

Z7T VFq 

where 1 denotes the direction which is perpendicular to q. One can see that the correction 
5 Im would be more singular than the free fermion contribution Im if q, O — > 
limit was taken with fixed Vt/vFq < 1. This result suggests that the perturbative expan- 
sion breaks down at low energies and the Fermi-liquid criterion are violated. Thus the 
gauge-dependent correction (which comes from the self-energy correction) to the trans- 
verse polarization function provides a similar picture as that from the singular one-particle 
Green's function [29]. 

Nevertheless, the perturbative corrections to the correlation functions should be cal- 
culated in a gauge-invariant way, thus one has to include the contributions from the vertex 
correction. The contribution to the transverse polarization function 5 Im 11^ coming from 
the vertex correction contains a singular term which exactly cancels the singular contri- 
bution from the self-energy correction. Thus the remnant terms in 6 Im 11^ provide the 
lowest order corrections to the transverse polarization function and have the following 
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form: 



6 Im II?! + 6 Im IT 



11 




7 



a 



(7^/x) 2/3 



+ 6 



(tQ/x) ' 

fcf,g 



(3) 



where a, 6 are dimensionless constants. One can see that the corrections calculated in 
a gauge-invariant way are always much less than the free fermion contribution as far as 
O <C vfQ and q <C k,F limits are concerned. Therefore, the perturbative expansion works 
quite well in this regime, at least up to the leading order gauge field corrections, and there is 
no need to go beyond the perturbation theory at this order. The validity of the perturbative 
expansion also indicates that the transverse polarization function is well described by the 
Fermi-liquid theory in the region of Vt <C vpq and q <C kp. This provides a very different 
picture from that obtained through the gauge-dependent one-particle Green's function. 

In this paper, we examine several gauge-invariant two-particle Green's functions or 
response functions in the limit of low frequency and long wavelength. It is shown that all 
the leading singular contributions from the self-energy correction are cancelled by the con- 
tributions from the vertex correction in systematic perturbation theories (which guarantee 
the gauge-invariance in each order of the perturbative expansion) . This cancellation is es- 
sentially due to the Ward identity. It is found that singular corrections to the two-particle 
Green's function do not appear for all ratios of Vt/vpq as far as the limit of low frequency 
and long wavelength limit is concerned. This kind of cancellation was also discussed by 
Ioffe and Kalmeyer [30] for a static gauge field. Recently, Khveshchenko and Stamp [23] 
performed non-perturbative calculations of one-particle and two-particle Green's functions 
using the so-called eikonal approximation. Even though they obtained a highly singular 
one-particle Green's function, the singularity does not show up in two-particle Green's 
functions for small q and O in this approximation. 

We also show explicitly that the gauge field propagator is not renormalized by the 
fluctuations beyond RPA up to two-loop order. Non-renormalization of the gauge field 
propagator was first discussed by Polchinski [28] in the framework of a self-consistent 
approach. In this approach, it is assumed that the dispersion relation of fermions is given 
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by u> oc ^ 2 (£k = k 2 /2m — and that of the gauge field is given by fi oc iq 3 , which are 
the results of one-loop corrections. Ignoring vertex correction by assuming the existence 
of a Migdal-type theorem, he showed that the assumed one-particle Green's function is 
self-consistent, and the polarization function is given by the same form as that of free 
fermions Im II^ = — (mv F /2n) (Q/vFq) for fi < ^/ 1 / 3 ^ 2 / 3 ^ 3 / 2 . As a result, the gauge field 
propagator is not renormalized because the dispersion relation of the gauge field is given 
by fi oc iq 3 . However, we would like to remark that his result is quite different from those 
obtained in this paper. One can check that the polarization function in the self-consistent 
approach has a different form compared to that of Fermi liquid for fi > ^/ 1 / 3 ^ 2 / 3 g 3 / 2 . 
However, in our perturbative calculation, the cancellation of anomalous terms from self- 
energy and vertex corrections leads to the result that the polarization functions have Fermi 
liquid forms for all q and fi as far as both are small. 

We have made several explicit calculations of two-particle Green's functions. In partic- 
ular, we consider a model given by Eq.(8) with u(q) = V /q 2 ~ v (v(r) oc V /r v , 1< rj < 2) 
which corresponds to the interaction between fermions in the problem of HFLL. We will 
present the non-analytic contributions (due to the gauge field fluctuations) to the two- 
particle Green's functions. The transverse polarization function iln(q, (]) up to two-loop 
order is found to have the following form. For fi <C vpq, we get 



lmn n (q,0) 



mvp fi 



2-7T VFq 



2 3 

mvp (t^Vx) 1 ^ , mvp (7^/x)~ 

J. (X ~ Fi 



7 



kp 



7 



k 2 F q 



(4) 



while for fi ^> VFq, 



Im n n (q,0) 



l + ?7 



3-V 



vf 7 1+n ^ 3 - r > 



8tt 2 (5 + 77) sin(^) rn x ih 



1 + c mvp 



3 / X N 1+ " 



7 



2^+3 

fi 1+1 



(5) 



where a, 6, c are positive dimensionless constants. 

The density-density correlation function n o(q, fi) is also calculated. We have a for- 
mula valid for any ratio of Q/vFq as long as fi and q are small (see Eq.(70)), but here we 
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just discuss limiting cases. For O <C vpq, we have 



im n 00 (q, fi) 



27T upg 



m O 



1 - 



4tt(5 + rf) cog ( r^i\ k F m x t$ 



1 + 1] 1 1 7 1 




On the other hand, for O ^> vfQ, 



Im n 00 (q, O) 



1 + 7? 1 




3-i? 

!+*7 3-7) 



) 



2 



(7) 



8tt 2 



(5 + ^)sin(^) 



Note that Im Ilnfg — > 0,f2) = -g— aim IIoo(g — ^ 0,O) is satisfied as it should be. Eqs.(4)- 

(7) are the main results of this paper. 

From the above gauge-invariant correlation functions, one can see that 
1) The corrections are irrelevant in the small q and Q limit regardless of the way how q 
and O approach to zero (for example, q — > limit may be taken first or O — > first, 
etc.). Therefore, non-perturbative calculations are not necessary. However, the sub- 
leading contributions are in general non-analytic due to the long range nature of the 
gauge interaction. The non-analytic sub-leading terms may have some experimental 
consequences. For example, the NMR relaxation rate 1/Ti in the problem of HTSC 
can be determined from n o(q, O). At low temperatures we have 



where n o plays the role of spin susceptibility in HTSC. Eq.(6) implies the following 
non-analytic correction to the free fermion result (only contributions from small q are 

1 5 + 17 

considered) ^ ^ ~ ^ T 1 + r i , where A is a constant and the first term is the 

result of Fermi liquid. Notice that this result is in disagreement with a result based 
on a renormalization group approach obtained in Ref.26, even near rj = 1. For HTSC 
i] = 2 and ^ 1 — ^. T 7 / 3 . Note that the non-analytic correction is very small so 
that the Fermi liquid form is preserved. 
2) q — > limit of the transverse polarization function indicates that the transport scat- 
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tering rate T tr (which determines the DC conductivity) scales as T tr oc 1 +-n at low 




q 
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frequencies (see Eq.(45) for more details). This result can be also obtained from the 

coefficient of the term which is proportional to q 2 in Im IIoo(q, O), and the relation 

o 2 

Im Tln(q — > 0,Q) = -^-Im U 00 (q — > 0,O). This result exactly agrees with those 
obtained by different approaches [12,16]. Note that r tr < O for 1 < 1] < 2. 

3) From Eq.(4), one can see that the gauge field corrections are smaller than the result 
of free fermions along the curve O oc q 1+r] which is the dispersion relation of the gauge 
field. Therefore, the gauge field propagator is not renormalized. As mentioned above, 
non-renormalization of the gauge field propagator was first discussed in Ref. [28] within 
a self-consistent argument. 

4) For i] < 2, the gauge field corrections to the polarization functions are less singular 
than the result of the free fermions for O < vpq. In particular, the edge of the 
particle-hole continuum in Im ITn and Im n o still occurs at O ~ vpq, where vp is 
finite and shifted from the bare fermi velocity as in the usual Fermi liquid theory. 
We conclude that the two-particle Green's functions are consistent with those of a 
Fermi-liquid with a finite effective mass. However, a combination of a divergent mass 
and divergent Fermi-liquid parameters cannot be ruled out. 

The remainder of the paper is organized as the following. In section II, we introduce 
the model and review some one-particle properties. In section III, the transverse polariza- 
tion function for q — > case is calculated. The cancellation of anomalous terms (coming 
from the self energy and the vertex correction) up to (l/iV) th order is explicitly shown 
(where N is the number of species of fermions). We also discuss the optical conductivity 
using the information of the calculated transverse polarization function. In section IV, we 
calculate the transverse polarization function for finite q <^kp case. It is also argued that 
the gauge field propagator is not renormalized up to two-loop order. In section V, the 
density-density correlation function is calculated up to two-loop order for finite q <C kp. 
In section VI, the results are compared to the conventional Fermi-liquid theory and their 
implication is discussed. We conclude this paper in section VII. 

8 



II. THE MODEL AND THE ONE-PARTICLE PROPERTIES 

The model is motivated by the above mentioned two strongly-correlated electronic sys- 
tems. It is constructed such that it includes the most important infrared singular behaviors 
of the one-particle Green's function. In this paper, we consider the zero temperature limit 
and use the Euclidean space formalism. The model in the Euclidean space is given by 



Z = J Dip Dip* Da^ e-f dT Sr C , 



where 



C = ip*(d - ia - n)ip - -^—ip*{di - iai) 2 ip + ia n f 

a 2 f ® 
+ y / dV(V x a(r)) v(r - r') (V x a(r / )) . 

Here v(q) = V /q 2 ~ v (v(r) oc Vo/r 77 , 1 < t] < 2), m is the bare mass of the fermion, 
and n/ is the average density of fermions. We choose the Coulomb gauge V • a = 0. 
Note that this model is incomplete for the problem of HFLL because of the absence of 
the Chern-Simons term. However, one may expect that it contains possible low energy 
singular behaviors because the most singular contribution to the one-particle properties 
comes from the transverse part of the gauge field. In the problem of HFLL, a = l/(27r0) 
and = 2 which is the number of flux quanta attached to the electron [1] . For the case of 
HTSC, one can take a = [12,13]. 

After integrating out fermions and including gauge field fluctuations up to one-loop 
order (RPA), the effective Lagrangian density of the gauge field is given by [1,12,13] 



--eff — - 



1 f d 2 q du> * . , . . . . . . 

2 J (2^)2 2^ a ^ q ' W ) D ^^ IUJ ) a M^) , (9) 



where 



D J~ ("a U° ll+ a 2 v(q)Q 2 ' ' (10) 



Here fx, v = 0, 1 and 1 represents the direction that is perpendicular to q. Uq and are 
given by the one- loop diagrams in Fig.l (a) and (b) respectively. In the limit of u> <^ vpq, 

9 



one can find that [1,12,13] 



2ir \ vfQ 

2n \u\ q 2 

— — H — 

kF q 24nm 



^ = ^ + ^— (11) 



a; 

7 — + XoQ 
Q 



2 



Therefore, the gauge field propagator can be expressed as 

D-i = -™ (l-M 
00 2tt V v F q 

D^ = ^ + ( X o + a 2 v(q))q 2 (12) 

M n 

~ 7 — + XQ V , 
Q 

where X = Xo + ct 2 Vo for rj = 2 and x = « 2 Vo for rj ^ 2. 

Since the longitudinal part of the gauge field is screened, the transverse part of the 
gauge field dominates the physics. The one-loop self energy correction due to the transverse 
part of the gauge field is calculated as (Fig. 2) [1,12,20] 



f d 2 q 


dv 


k x q 


1 (27T) 2 


2^ 


m 



2 



E(k, iu) = 

~ — i A Icul 1 ^ sgn(o;) , 

where 

A = 



G (k + q,zw + iv) Dn(q,zi/) 



47r sin(^p-) 7^+! x 1+?7 



and 1 (k, zcu) = zcu — £k (£k = ^ — A 4 )- The self energy as a function of real frequency O 
(in the Minkowski space) can be obtained from the analytic continuation of S(k, ioj), i.e., 
S(k, O) = E(k, iuj -> O + z<5). Note that |Im E(k, 0)| oc > |f2j for sufficiently small 

O or |0| <C A^ 1 . Therefore, the quasi-particle (the dressed fermion) is not well defined. 

This can be also seen from the spectral function of fermions. The spectral function 
can be obtained from the imaginary part of the retarded Green's function: A(k, O) = 

Im G R (k, O) = -1 Im G(k, iu ^Q + i5), where G'^kju) = G^Ql, iw) - E(k,iw). 
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In the low frequency limit, 



A(k, n) « i ^ " " , , (14) 

7T ( Al _^)2 + ( A2 — )2 



where Ai = A cos 



7T / y-1 
2 U+l 



7T ( 77- A 

2 



Note that the maximum of 



and A2 = A sin 

A(k, 0) appears at O ~ ( ^ J . However, the width of the broad peak is also order AO ~ 

, , its 

( • Therefore, the Landau criterion for the existence of quasi-particles (AO <C O) 
is marginally violated. 

If we assumed that there is a well-defined Fermi wave vector kp = (4-7rn f) 1 ^ 2 and tried 
to fit the result to the usual quasi-particle picture, the energy spectrum of the quasi-particle 
would be [1] 

I+T7 

ek oc \k — kF\ 2 (15) 
for k sufficiently close to kF- From ^ = %r\ k=kp , the effective mass diverges as 

m oc \k — kF\ 2 oc |ek| 11+1 ■ (16) 

This suggests that at least some modifications to the conventional Fermi-liquid theory are 
necessary as far as the one-particle Green's function is concerned. 

There have been also some non-perturbative calculations of the one-particle Green's 
function [21-24], which were motivated by the singular perturbative correction at low en- 
ergies. The results look very different from that obtained by the lowest order perturbative 
calculation and even exponentially decaying one-particle Green's function is found in the 
so-called eikonal limit [23]. 

From these results, one may doubt the validity of the quasi-particle picture although 
a modified Fermi liquid description is proposed for the case of the HFLL [1]. However, 
one should also remember that the one-particle Green's function is not gauge invariant. 
This can be easily seen in the path integral representation of the one-particle Green's 
function [12,21] of a fermion interacting with a gauge field, i.e., each path acquires a phase 
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factor e §o dt &{v,t ) dr/dt ^Yada. is manifestly not gauge invariant. Therefore, it is very 
important to examine gauge-invariant quantities. As the first example, we will calculate 
the polarization function for q — > case in the next section. 



III. THE TRANSVERSE POLARIZATION FUNCTION FOR q -> 
AND OPTICAL CONDUCTIVITY 

Let us consider a large N generalized model of Eq.(8), where N is the number of 
species of fermions. In this model, each fermion bubble carries a factor of N and each 
gauge field line gives a factor of 1/N. Thus, for example, IIq and IL^ obtained in the 
previous section should be multiplied by N. 

In this section, we consider only the q — > case of the transverse polarization function: 
nn(q — > 0, if). However, the relevant diagrams are the same even for q ^ case. The 
leading order contribution is which is proportional to N. The relevant diagrams in the 
next order (i.e. (l/AT)°th order) are given by Fig. 3 (a)-(g). For convenience let us define 
the following quantities: = (a) + (b) and II^ = (c) + (d). The formal expressions of 
these quantities for q — > case are given by 

d 2 k doo TA; 2 -(k-q) 2 ^ 



n 



(i) 
ii 



(2tt) 2 2tt 



S(k, ioo) [G (k, iu)Y G (k, ioo + iv) , (17) 



and 



n 



(2) 
11 



d 2 k doo 



k 2 - (k • q) 

m 2 



E(k, ioo + iv) [G (k, ioo + iv)] 2 G (k, ioo) . (18) 



(2tt) 2 2tt 

These two equations can be rewritten as 

d 2 k doo fA; 2 -(k-q) 2 l E(k,ioo) 



n 



(i) 

; : ' (2tt) 2 2tt [ %v 

x ([G (k,iw)] 2 - G (k,ioo) G (k, ioo + iv )) , 

2 



(19a) 



(2) _ /" crfc cL; 

Ul1 " 7 (2tt) 2 27T 



k 2 - (k • q) 2 



E(k, icu + if) 



iv 



(19b) 



x ([(?o(k, iw + if)] 2 - G (k, icu + iv) G (k, ioo)) 
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If we add (19a) and (19b), the first terms in each polarization bubble are cancelled by each 
other and the remaning parts give us 

d 2 k du Tfc 2 -(k-q) 2 ' 



n (i) , n (2) _ 
n n -t- — 



(2tt) 2 2tt 



E(k, iu) - E(k, iu + iv) . s 

x — — - G (k,iu) G Ck, iuj + iv) 

iv 



(20) 



From the above expression, it can be easily seen that the contributions from (b) and (d) 
are automatically cancelled because the self energy corrections in these diagrams are just 
the same constants. 

Next we consider the diagram given in Fig. 3 (e). Here we have to include the vertex 
correction for q — > case (Fig. 4): 

;a2" 



ri(k, q — > 0; iu, iv) = 



f d 2 q' 


dv' 1 


' k 1 + q[\ 


k 2 - (k • 


1 (27T) 2 


2^ \ 


m ) 


m 2 



(21) 



x G (k + q', iiv + iv') G (k + q', iu + iv' + iv) Dn(q', iv' 



Then nff(q^ 0, zV) can be written as 



n (3) _ 



r d 2 k du 


'_ki 


J (2tt) 2 2tt 


m 



Ti(k, q — > 0; iu, iv) Go(k, iu) Go(k, iu + iv) 



— n 11 -t- n n , 



(22) 



where 



ii n — 



d k du 



x 



(2tt) 2 2tt 
d 2 q' dv' 



(2tt) 2 2tt 



fc 2 -(k-q) 

m 2 

k 2 - (k ■ q') 2 

m 2 



Go(k, iu) Go(k, iu + iv) 



(23) 



and 



x G (k + q', iu + iv') G (k + q', iu + iv' + iv) Du(q' , iv') 
d 2 k du 



n^' 2) = - J ^ G (k, iu/) G (k, iw + zV) 



x 



dV dv' (q^k x 



) 2 2tt 



k 2 - (k-q') 2 



(24) 



x G (k + q'jicj + zV) G (k + q', iu + iv' + iv) Du(q' , iv') . 



(3 1) (3 2) 

Here we would like to point out that Tl\± is more singular than . This can be 



easily seen from the fact that n^' 2 ^ can be obtained by replacing k\/m 2 = 

13 



fc 2 ~(k-q) 2 



in the integrand of Eq.(23) by q[ki/m 2 . Using q[ = sin 6>k q + q± cos 6>k q and £k+q ~ 
£k + vfQ\\ + q±/2 m i one can do the integrals over q'^ and q' ± in Eq.(24). Since the con- 
tribution from q'j_ cos 6>k q term becomes an odd function of q' ± , this term vanishes. By a 
formal manipulation, one can replace q'^ by q'±/kF so that q[ factor becomes effectively 
(q' 2 ±/kF) sin 6>k q - Since the integrand is dominated by \v\ ~ {x/l) scaling given 

by the pole of the gauge field propagator, replacing k± by q[ gives rise to an additional 
factor which is proportional to |z/|i+^. Therefore, n^' 2 ^ should be less singular than n^' 1 ^ 

2 

by the factor \v\ 1 + t » in the low frequency limit. 



(3 1) 

Note that 11-^' can be rewritten as 
n (3,i) _ _ / d2k 

11 ~ J (2tt) 2 27T 



k 2 - (k • q) : 



To(k, q — > 0; iu, iv)Go(k, iui) Go(k, iiv + iv) , 

(25) 



where r is the scalar vertex: 

r (k, q; iu.iv) = 

> 

From the relation, 



f d 2 q' dv' 


~k 2 - 


-(k-q') 2 ~ 


I (2tt) 2 2tt 




m 2 



(26) 



x G (k + q', iuj + iv') G (k + q' + q, ioo + iv' + iv) Dn(q', iv') 



v,/, • s v,/, • ■ s f d 2 q' dv' 
E(k,z W )-E(k,z W + zz,)=y ^2 2^ 



fc — (k • q') 



>\2 



x [G (k+q',iu; + ii/) - G (k + q' , iu + iv' + iv)] £> n (q',iz/) 



I 



d 2 q' dv' 



(2tt) 2 2tt 



A; 2 - (k • q') 2 



iv 



x G (k + q',zo; + iv') G (k + q', ioo + iv' + iv) Dn(q', zz/) . 



we get the following identity: 

E(k, iu) — E(k, zcu + zV) 



r (k, q — > 0; iuj, iv) 



(27) 



(28) 



This is nothing but the Ward identity. From Eqs.(20), (25), and (28), we have 

nS?+nS?+nSf = 0. 



(29) 



14 



(3 2) 

Now the remaining piece is just 11^' . Following the procedures of integration mentioned 
above, in the low frequency limit, we get 

1 + 77 



n (3,2) 
Ai ll 



3-V 

vp 7 1+r < 



4tt 2 (5 + rj) sin I f 



n 



3-V 

\v\ 1+r > 



TTl l + r) 



(30) 



Here it is worthwhile to compare this result with Tl\^ +11^ and IT^' 1 ^, i.e., the results 



before cancellation. By a straightforward calculation, one can get 



n (l) , tt(2) 



(3 1) 

In order to calculate H\{ , the vertex correction should be calculated. The vertex correc- 



2 (l + r/) 
7T (3 + rf) 



77-I 



m v P A \v\ r i+ 1 



(31) 



tion r (k, q — > 0; iu;, zV) is found to be 



1 



1 



7 2. sin (^) - 



\ UJ \l\ 1+71 



X 



sgn(w) - 



+ ^ It \ 1+ " 



x 



sgn(cu + i/) 



(32) 



f3 1) 

Using Eqs.(25) and (32), can be calculated as 



n (3,l) 
Ai ll 



m v. 



l + r? 

2^sm^) V3+^ 



2 

^77 + 1 1+77 



77-1 

\v\ T >+ 1 . 



(33) 



(1) (2) (3 1) (3 2) 

Note that, as mentioned above, II-q + II-q and 11^' are more singular than H\{ by 

_ 2 

\v\ 1 + r > in the low frequency limit. The important point is that these singular terms are 
cancelled by each other due to the Ward identity. 

Now let us look at the diagrams of (f) and (g). Let = (f) and II ^ = (g). The 
formal expressions of these diagrams for q — > case are given by 

j2„/ a,, i a,. j fly duo" 



(4 ) _ f d 2 q' c/V d 2 k' du' d 2 k 
11 ~ J (2+p 2^ (2^)2 "27 (2^ 



x 



k' ■ k" - (k 7 ■ q') (k" • q') 



rrr 



) 2 2tt 

Dn(q',ii/') Dn(q',ii/' + ii/) 



(34) 



x G (k / , iu/) G (k', iu/ + ii/) G (k' - q', iu/ - iv') 
x G (k", iu") G (k", iu;" + iv) G (k" - q', iu" - iv') , 
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and 



n (5) _ 



X 



d 2 q' du' d 2 k' du' d 2 k" du" 



(2tt) 2 2tt (2tt) 2 2tt (2tt) 2 2tt 



k' ■ k" - (k' ■ q') (k" • q') 



Dn(q',ii/') Dn(q',M/ + iz/) 



(35) 



x G (k', zu/) G (k', z'u/ + zV) G (k' - q', z'u/ - ii/') 

x Go(k", iu") Go(k" , iu" + iu) Go(k" + q', iu" + iv' + iu) . 
By changing variables as q' — > — q', v' — > —v' — u and using Dn(-q', — ii/) = L> 11 (q / , zV), 
we get 

72 V 



x 



<i 2 /c du ki f k sin #k q < x 



(27r) 2 27r m \ m 



G (k, iu) Go(k, iu + iu) 

2 



(36) 



x (Go(k + q', iu + iv' + iu) + G (k - q', iu - iu')) 



where 6>k q ' is the angle between k and q'. In the low frequency limit, we get 



TTl 1+17 



(37) 



where c\ is a constant. One can also show that 



n 



m v. 



(4) 

11 c i — -g— 



3-1) 



11 ~ c 



m v, 



V-l VF 7 1 + ,) 3-7? 



(38) 



T?-l 2 

'y 77+1 1 + r ' 



f 1+ 1 — C\ — 



TO 1+77 



where c is a constant. That is, there is also a cancellation between the singular parts of 
Gathering all the previous informations and using n^q — > 0, zV) = we can 

(39) 



nff and nff 



conclude that 



3-77 



n 



Nn kp 7 1 + ? ' 3-77 



ii 



c 2 



9 4 
777/ ^T+^ 



Z7 



up to (l/AT)°th order, where C2 is a constant. 
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In order to calculate the optical conductivity, we have to consider the bubble diagrams 
with two external lines that represent the coupling to the external vector potential while 
the internal gauge field lines are due to a M . There are additional diagrams generated by 
ijj^ a^A^ip vertex. All the additional diagrams except one (shown in Fig. 5 (a)) vanish due 
to the symmetry of the integrand. A typical diagram which vanishes is shown in Fig. 5 (b). 
It turns out that the diagram represented by Fig. 5 (a) gives an imaginary part which is 
higher order in frequency compared to |z/|i+^ so that it is irrelevant in the low freqency 
limit. Now we can use the imaginary part of the transverse polarization function in the 
Minkowski space ITii(q — > 0, O) = ITii(q — > 0, iv — > O + iS) to calculate the real part of 
the optical conductivity: 



2 Im n n (0) 



Re <t(Q) = -e 2 . (40) 



From Eq.(39), Re (x(fi) is given by 



,2, 5=2 



Re a(O) oc ^ ^ 0" 2 (^) 



, 2 s - ■ (41) 

X 1 



If there were no cancellation, the result would look quite different. For example, if we did 
not consider the vertex correction, the result from II^ + 11^ would be 



Re a nv (Q) oc ^ ^_ O"^ , (42) 

where a nv represents the conductivity without vertex correction. 

Now we are going to show that the right answer given by Eq.(41) is consistent with a 
modified Drude formula if we assume that the transport scattering rate (which is the inverse 

1 1 3 — ?7 4 4 

of the transport time Tt r ) of the fermion is given by r tr (0) oc (7!+^ /x 1+v ) Q 1+ri ■ 

First of all, for later convenience, let us calculate the inverse of the transport time r t ° r 
of the fermion [12] using the imaginary part of the self energy E(k, O). For this purpose, 
we can just include the factor 1 — cos = 2 sin 2 (0/2) in the integrand of the expression 
for Im E(k, fi), where is the angle between the wave vector of the fermion and that of 
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1 

the gauge field [12]. Using the fact that sin (6/2) ps q/2kp and q ~ ^ ms ide the 

integral [12], we get 

3-V 

I I I 71+1 4 

— « ^7 T ^ fiTT " ( 43 ) 
r t u r iv m/Cf 

Therefore, we will essentially show that our result of the optical conductivity is consistent 
with the identification of I\ r = 1/ r t ° r or r tr = T t ° r in a modified Drude formula. 

The Drude formula that is appropriate to the large N generalized model is given by 

In the large iV limit, if we assume r tr = l/r t ° r oc 1/N, 

Re a O « -£ oc -i-j- O 2 ^ J . (45) 

This is the same result as that of Eq.(41). The result of Eq.(42) can be reproduced in the 

. 77— 1 _ 2 2 

same way if we assume that r tr (0) oc {mv F ) (7 '?+ T x ~ ) which is essentially the 
imaginary part of the self energy E(k, O). Therefore, the optical conductivity is consistent 
with the choice of l/r t ° r rather than just the naive scattering rate (given by the self energy) 
as the transport scattering rate. Since the singular contribution, which gives Eq.(42), is 
cancelled by the vertex correction, we can again say that the leading singular behaviors of 
one-particle properties do not show up in the optical conductivity. 

For finite temperature, one can replace O by T in r tr . Note that the DC-limit of the 
optical conductivity Re cr(0 — > 0) = Nr //f jj- cannot be obtained by the 1/N expansion. 
However, one can infer the DC-limit by assuming that the full Re cr(O) is given by Eq.(44) 

(with r tr = r tr (T)) which is consistent with the result of the large- N limit of the optical 

4 4 
conductivity. If r tr oc T'+i was used, one would get Re cr(T) oc T ^ [12]. One the other 

2 2 

hand, one would get Re a nv (T) oc T i+v if r tr oc T'+i was used. In Ref.[19], the authors 
concluded that the resistivity of the system is proportional to T 2 / 3 for the short-range 
interaction (77 = 2) and this is consistent with the latter case. Therefore, our result is in 
disagreement with their conclusion about the resistivity. 
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IV. THE TRANSVERSE POLARIZATION FUNCTION FOR FINITE q < k F 
AND NON-RENORMALIZATION OF THE GAUGE FIELD PROPAGATOR 

It is not easy to find the polarization function for arbitrary q and v. However, some 
simplifications can be made for q <C k F case. In this section, we calculate IIii(q, iv) for 
finite 5 < fcf up to two-loop order. We set N = 1 first, and discuss the extension to the 
large- N case later. 

First of all, II ^ and for finite q have the following form: 



(i) _ _ f d 2 k du 
11 ~ J (2-) 



r) 2 2tt 



n ( 2 )__ [*> 

11 - J 



d 2 k du 



) 2 2tt 



A; 2 -(k-q) 2 

m 2 

k 2 - (k • q) 2 

m 2 



E(k, iu) [G (k, iu)} 2 G (k + q, iu + iv) , 



E(k + q, iu + iv) 



(46) 



x [Go(k + q, iuj + iv)] 2 G (k, iu) . 
Using the similar method as that used in section III, one can obtain 



n (l) , tt(2) 
iA ll + 



d 2 k du 



(2tt) 2 2tt 



(k-q) 



m 1 - 



G (k, iu) G (k + q, iu + iv) 



E(k, iu) — E(k + q, iu + iv) 



(47) 



iv — v F q cos 6*k q 
Next we should consider the vertex correction (Fig. 4) for finite q: 

r d 2 q' dv' 

ri(k,q;ia;,ii/) = J — A(k,q,q') B(k,q,q') 



(48) 



where 



A = 



x G (k + q',iu + iv') G (k + q' + q, iu + iv' + iv) -Dii(q', iv') , 
k 1 + q' 1 +q 1 /2 k x + q[ 



m 



m 



1 



B = — 

m 



(k + q'/2) • (k + q + q'/2) - (k + q'/2) • q' (k + q + q'/2) • q' 
For q <C k F and |k| pa k F , the following approximation can be made 

A; 2 -(k-q') 2 



(49) 



B 



m 



(50) 



Using this approximation, one can show that 

tt(3) f d 2 k du ki~\ . . , . . , 

llii = — / , .„ — 1 i(k, q; iu, iv) G (k, iu) G (k + q, iu + %v) 



n (3,3) n (3,4) 



(51) 
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where 



n (3,3) / d 2 k du 

11 " J (In)* 2n 



^-(k-q) 



To(k, q; iu, iu) Go(k, iu) Go(k + q, iu + iu) , 



n ( 3 > 4 ) - 



x 



d 2 k du 



(2tt) 2 2tt 



Go(k, iu) Go(k + q, iu + iu) 



d 2 q' du' ( q[k 



(2tt) 2 2tx \m 



k 2 — (k • q 



A2 



x G (k + q', iu + iu') G (k + q' + q, iu + iu' + iu) Dn(q', iu') . 



(52) 



First, let us calculate the scalar vertex part r (k, q; iu, iu). We use £k+q' ~ £k + 
VFq\\ + q'±/2m and 6c+ q '+ q ~ £k + v F q'^ + v F q cos 6> kq + ^ sin 6> kq + q' 2 L /2m (where 
<?ll = q' cos ^kq' and q' ± = q' sin 6>k q ') to perform the integral in Eq.(26). Using the 
fact that the important region of q' is the order of u~ C 1 so that q'/k m q'/k F <C 1, 
we conclude [23,27,28] that q'^/k F (q' ± /k F ) 2 and we can approximate the gauge field 
propagator as Dn(q',w') ps l/(7|z/|/|<j^_| + X I 1 77 ) - After performing q'^ integral, we get 

T (k, q; iu, iu) ss -zv F j ^- J ^ (sgn(w + u') - sgn(w + v + i/)) 



(53) 



x 



iu - v F q cos # kq - ^ sin # kq 7^j + xW±\ v 



Now z/ integral gives 



T (k, q; iu,iu) 



7 J-fc F 



_L 



X 



In 1 + 



fc_F f + u>i?g cos 6>k q + i-^ sin 6>k q 

sgn(w) - In ( 1+ ,^ n-TJ 7 . I + l/ ) 



By changing variables, one can get the following formula. 

T (k,q; iu,iu) 
v F 1 



(54) 



1 



7 TV 2 U -\- W F q COS $kq 



X 



u 
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X 



1 + 7J 



(g/m) sin 6> kq 

I^Q COS 6*kq — if 



H7 



|w + z/|7y+" / (q/m) sin # kq 

.r \ u -\- u, 



X 



v F q cos 6*k q — if 



sgn(w) 

|cu + z/|7 

X 



(55) 



1+7, 



sgn(u; + u) 
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Here F(u,x) is defined as 



F{u),x) = I dy \y\ - 
•J -y c 



1 + xy 



(56) 



i 

1+7, 



where y c = (^]^J • ^ can t> e easily shown that q — > limit of Eq.(55) is given by 
Eq.(32). On the other hand, the self energy can be rewritten as 



£(k,u;) 



_.Vf_ (lulY} 1+ " 



7T 2 7 



X 



(57) 



Collecting these results, it can be shown that 



n£> + n£) + nfr 3 > 



<i 2 /c du; 



x 



(2tt) 2 2tt 
w.f 1 



A; 2 - (k • q) 
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G (k, iu;) G (k + q, iu + iv) 



TV 2 J VFq COS 6>kq — iv 



I(u) - I (a + v) 



(58) 



where 



X 

F l w 



1+7, 



sgn(w) 
(<?/m) sin 6> kq 



X 



1 + 7? 



F(w,0) 



(59) 



vpq cos 6*k q — if 

The integrals in Eq.(58) can be evaluated as the following. Using J d 2 k/(27r) 2 = 
(m/27r) / d£k / d6>kq/27r, one can perform £k integral easily. The angular integral over 
6*k q can be done by contour integration, which requires long algebraic manipulations. The 
remaining u integral and the y integral in I (to) of Eq.(59) can be evaluated by scaling 
the integration variables and expanding the integrand in some limits. More details of the 
calculation will be demonstrated in the later evalution of the density-density correlation 
function (see the discussions about Eqs.(68)-(70) in section V) which can be more easily 
calculated. First, for \u\ <C VFq, 



nff+ng' + ng^ca 



tu Vp \u\ (l\v\/x) 1+r > 



7 v Fq 



fee 



(60) 
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while, in the other limit \v\ ^> vpq, we get 



7 M 



(7/x) 1+ " 



777, 1/ 77+1 



(61) 



where C3 and C4 are dimensionless constants. 

The calculation of ng' 4 ' ) can be also done by the similar method used in the evaluation 



of ng' 3 ' ) . First, for \u\ <C vf<1, we get 
r(3,4) , 



n 



11 



m Vp \v\ 
7 v Fq 



(7lH/x) 1+ " , hW\/x) 1+ " 

C5 To h C6 



(62) 



whereas, in the other limit ^ vpq, 



n 



(3,4) 
11 



1+77 



1 



3-1? 

t) F ^1+17 3-7) 

— — W\ 1+r > - c 7 



2„,3 
F 



v F q 



7 fn 2 (x/l) 1+v M 1+v 



, (63) 



where C5, C6 and C7 are dimensionless constants. 



(1) (2) (3 3) (3 4) 
From the above results, it can be shown that in^+n^+n^l < \U^{ > \ for relevant 

limits. Therefore, the imaginary part of the transverse polarization function IIii(q, O) (in 

the Minkowski space) up to two-loop order is given by the following formulae. For O <C vpq, 

we get 

r 2 3 



Im n n (q, O) 



mvp O 
2tv vpq 



1-a 



7 



7 



kin 



(64) 



where a and b are dimensionless constants. Note that the correction is small as far as 
1 < V < 2 is concerned. On the other hand, for O 3> vpq, we have 

I + 77 1 



3-rj 



Im n n (q,0) 



Dp O/l + rj 3-7, 



87r 2 (5 + ?7) sin m x tt^ 



-1 3 I X 

1 + c mvp 

' ' 7 



1+77 



22+3 



(65) 



where c is a dimensionless constant. 

For O > vpQ, there is no contribution to Im Iln from the free fermion bubble because 
the regime is outside the particle-hole continuum. Therefore, any non-zero contribution 

22 



to Im Iln for J] » « F g entirely comes from the gauge field correction. Note that the first 

l+?7 1 2t; + 2 

term in Eq.(65) dominates for Q > (mvp) 2r >+ 3 (x/l) 2r >+ 3 q 2r >+ 3 . On the other hand, the 

1+77 1 277+2 

second term becomes more important for vfQ <C O < {mvp) 2r >+ 3 (x/l) 2v+3 q 2r > +3 so that 

2 -Tl 

A I + 77 2 

Im iln oc vp— 3 in this regime. As we approach the line given by Q = vfQ, Im Iln 

X 1 + r > fiTT^ 



4+2 2-2 , , 



becomes vp +v 7 ^ q 1 +^ as a function of q. 



2 

In the case of O <C i>f</S the free fermion bubble gives Imn n = Note 



that Im IIii(q, O) 



27T vpq 



<> 

27r t> F q 



1 _ m« F 



2-k vp-q' 

2 " 



3 

mv F (TjVx) 1+T? 
7 fcjq 



for n < (x/7)q 1+?? and Im n n « 
for (x/l)q 1+v < ^ ^ vfQ- It is gratifying to note that, 



along the line Q = vpq, the correction to Im Iln given by the above expression agrees 
in its q dependence with that obtained by approaching from Q ^> vpq given in the last 
paragraph. In any case, the corrections are small compared to the free fermion result for 
1< TJ < 2. 

Using the result of iln for \u\ <C vpq, we can discuss the issue of the renormalization of 
the gauge field propagator. Recall that the dispersion relation of the gauge field obtained 
from the one-loop correction is given by \u\ ~ (x/l)Q v [1,12,13] , which is below the 
line of \v\ = VFq for sufficiently small q. Along the line of \u\ ~ (x/l)Q 1+r> ? one can 
easily see that the correction to II n is smaller by In ^ £ - (j^ j . Therefore, the gauge field 
propagator is not renormalized up to two-loop order. As mentioned in the introduction, 
non-renormalization of the gauge field propagator was first discussed by Polchinski within a 
self-consistent argument and without vertex correction. In Ref.[19], the authors discussed 
the relevance of Y^\a^) and r^(a M ), which are coefficients of the a 3 and a 4 terms in 
the expansion of the effective action of the gauge field. They concluded that T^\a^) 
and r( 4 )(a M ) are irrelevant so that the gauge field is not renormalized. Since the two-loop 
diagrams we considered are generated from r^(a M ), our calculation is consistent with their 
conclusion. By analogy, we expect that II ^ and U^±± are irrelevant for the renormalization 
of the gauge field because these are generated from Y^\a^). We also directly evaluated 
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r( 3 )(a M ) and confirmed the argument of Ref.[19]. Therefore, one can expect that the gauge 
field is not renormalized up to (l/iV)°th order in the 1/N expansion. That is, the RPA 
calculation gives the leading contributions in the low energy limit. 

V. THE DENSITY-DENSITY CORRELATION FUNCTION FOR FINITE q < k F 

The polarization function for the density channel n o(q, can be also calculated in 

a similar way as used in section IV. In this section, we consider the two-loop corrections 

given by Fig. 3 (a)-(e) and finite q <C k F case. The sum of the contributions from the 

self-energy corrections given by Fig. 3 (a)-(d) can be written as 

fu f d 2 k dui _ . . „ . . . E(k, iu) — E(k + q, iu + iv) 
n oo » / 77) \2 V~ Go k, iu) G k + q, iu + iv) 1 / ^ \ , 66 

J (2/Tj 2 2% IV — Vpq COS tfkq 

while the contribution given by Fig. 3 (e), which comes from the vertex correction, can be 
also written as 

n oo = ~ J (7^2 ^ r °( k ' ^ iuJ > iv ) G o( k ' iu) ) G o(k + q, iu + iv) . (67) 

Using Eqs.(55) and (57), it can be shown that 

n$ + n£ } « - J -J^ ^ G (k, iw) G (k + q, + «/) 

i 

J(a;) - J(u; + i/) 



X 



(68) 



7T 2 7 COS 6>kq — ^ 

where /(u;) is given by Eq.(59). Using J d 2 h/{2n) 2 = (m/2n) J <i£k J <i6>k q /27r, one can 
easily perform £k integral, which generates the additional factor vpq cos #k q — iv in the 
denomenator of the integrand of Eq.(68). Recalling that I{u) also has an angle dependence 
6*k q , one can perform the angular integral over 6>k q by contour integration, which requires 
long algebraic manipulations. After rescaling the u integral by a new variable x and the y 
integral in I(u) (see Eqs.(56) and (59)) by newly defined y, we get 



tt(d t.(2) 24 \v\ r , r , , / a/?^ 



7r j 7^r jo jo v y 

(69) 

ctl la] 



x 



(l + a 2 )^/l + a 2 +y 2 (1 + a 2 ) 3 / 5 
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1 

where a = and 3 = ( ) 1+v . 
integrals can be done, yielding 

ai \a\ 



In the small frequency v limit, the parameter 



n (1) + n (2) - 

il 00 + il r 



l 00 



ky\ (1 + a 2 ) 3 / 2 

l + rj 1 11 

4tt 2 (5 + rf) gin ( kp'j v F q 



7rl 
X 



4 

1+77 



or 



(70) 



(l + CK 2)5/2 ' 



\1+V J 

where a\ is an undetermined constant. This formula is valid for all ratios of q and z/, as 
long as both are small. Note that the first term gives only an analytic contribution, which 
also arises in the usual Fermi liquid theory. Similar methods can be used to produce a 
somewhat more complicated formula valid for all a for the transverse polarization function 
rin (for example, Eqs.(52) and (58) can be evaluated by a similar method). 

After dropping the analytic contribution, we combine the free fermion contribution 
and perform analytic continuation to get, for O <C i>i?g, 



imn o(<?, fi) « -7^- 



2% vpq 
and for O ^> VFq, 

im n 00 (<?, O) 



1 



1 + 77 



1 



4tt(5 + 77) CQS (^tt) X 1 ^ 



3-77 . 
1 7 / O \ 
' — — O.1+77 



v F q 



J 



(71) 



I + 57 



3-r, 



8tT 2 (5 + Tj) gin ( _2_n_^ /C F X T+^ 



1+t? 



1 ryl + r, 3-77 fVpqY 



(72) 



Note that Im Tln(q — > 0, 0) = -|^Im n o((Z — > 0,O) is satisfied. Therefore, both of 
Im Hn(q — > 0, O) and Im IIoo((Z — *• 0, O) give the same answer for the optical conductivity 
given by Eq.(41). 

VI. COMPARISION TO THE FERMI LIQUID THEORY 

In section III, it was shown that the resulting conductivity is consistent with a mod- 
elled Drude formula. In this section, we try to fit this result to the Fermi liquid theory 
framework to extract informations about the Fermi liquid parameters and examine whether 
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the gauge field induces some singular or divergent parameters. In the Fermi liquid theory, 
the conductivity for N species of fermions is given by [31] 

Nne 2 r 

= — i 7w — / — *Y ' ? 3 

m* 1 — iilr(m/m*) 

or 

Nne 2 Ft 

Re '<"> ~ ^ ■ (74) 

where r tr = Tsc 2 ^-, Tsc = 1/r is the scattering rate and r is the scattering time. Here m* 
is the effective mass of the fermion. Using the fact r tr oc 1/N in the large N limit, we get 

Comparing the above result with Eq.(41) which is a result of the 1/N expansion, we can 
again identify T tr with l/r t ° r given in Eq.(43). Therefore, we can conclude that r tr = Tsc^- 

4 

scales as after including 1/N corrections due to the gauge field fluctuations. 

In the following we will directly compare our perturbative result for n o with the 
density-density correlation function in the Fermi liquid theory. Our goal is to find out 
whether the perturbative result can be consistent with a Fermi liquid theory made up of 
quasi-particles with a divergent effective mass m* as suggested, for example, by Eq.(16). 
First we consider the limit O = 0, q — > 0, where it is well known that the Fermi liquid 
theory predicts 

n in :»!!-»)- n; ( q - 0. j] = 0) 

n °" (q 0,Q - 0) -i + /o,n; o(q ^o,si = o)' (76) 



n°-n° 



where IIq = — f p—ji o-f e * -e** ) * s ^ ne ^ ree f erm i° n response fuunction with an effective 
mass m* and /o s is the angular average of the Fermi liquid interaction parameter f pp > . In 
two dimensions, for small q limit, 

n; o(q , „) = _;!(!__£= ^ - 1) + o(i - *>)) , ,77) 

where x = VL/v* F q. In Euclidean space, the above formula can be reduced to 

m* / \a\ 



n So(q >) = --^--^= ) , (78) 
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where a = vjv F q. Since IIo (q — > 0, O = 0) oc m*, the fact that n o(q — > 0, O = 0) 
is not enhanced implies that /o s is a finite constant. However, this does not imply that 
the leading order term in the perturbative expansion of fo s is finite. In fact, it is clear 
from an expansion of Eq.(76) that if the leading order correction to m is singular, then 
the contribution to /o s at the same order should be also singular since n o has no singular 
correction in the lowest order perturbation theory. 

Next we consider the full q, O dependence of n o for small q and O. We are motivated 
by the belief that, in the Fermi liquid theory, Im n o(q, should exhibit the edge of 
the particle-hole continuum along the line O = v F q. However, when 11 ^ 0, a simple 
formula such as Eq.(76) does not exist for n o(q, O). In particular, n o(q, O) in general 
depends on the higher moment angular average of the Landau functions, and not just /o s . 
Nevertheless, the Fermi liquid theory makes a precise prediction for n o(q, for all q, O in 
terms of m* and the interaction parameter f pp i . This is given by the quantum Boltzmann 
equation for the quasi-particle distribution function n p = n p + 5n p in the Fermi liquid 
theory, where n p is the distribution function for the free fermion system with an effective 
mass m*: 



O- (e 



p+q/2 fc p-q/2 



5n r 



~ ( n p+q/2 - n p-q/2/ 



d 2 p' 



/p P ' <5n p /(q,fi) 



= . 



(79) 



Here e* is the quasi-particle energy, t/(q, O) is the external potential, and f pp > is the 
Fermi-liquid interaction parameter. The linear response of 5n p to the external potential 
can be calculated from Eq.(79) (to the first order in / pp ') : 



Sn p (q, O) 



r d 2 P ' ,. 

p J (27r) 2 tpIpp " p 





(2tt) 2 " pJpp p 



U(q,n) 



p+q/2 p-q/2 



(80) 
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p+q/2 " 
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e P -q/2) 



,2 

The change in the density of the fermions <5p(q, O) = J j2~ji <5™p(q 5 ^) is given by 
5p(q,H) 



<i 2 p n p - n p /" cFp cFp 



p -n p - q , /^V f (81) 
(2tt) 2 n - (e; - e p _ q ) + 7 (2tt) 4 C p^' C p' + ' " " ' 

where • • • represents the higher order terms in f pp > . The second term is just the diagram 
given in Fig. 3 (e), but with a frequency independent interaction / pp '. 

Let us now examine what happens to the edge in the particle-hole continuum according 
to our perturbative results. The gauge interaction may induce non-zero Fermi-liquid inter- 
action function f pp > and a change in the Fermi velocity 5vp. From Eq.(78) and Eq.(81), 
a change in the Fermi velocity 5vp and the appearance of the Fermi liquid interaction 
parameter induce the following change in the density-density correlation function: 

ATT 5VF ( TT* kF ^ \ f d2p ' n f n 

6U °° = - V r n °° + 2^(1 + ^)3/2; " J ~^yr Cp/pp,Cp ' • (82) 

If we assume a power law behavior for f pp > ~ [p^ijvp with A < 1 (i.e., finite /o s ), one can 
show that the second term in Eq.(82) cannot produce the singular term (1 + a 2 ) -3 / 2 near 
a 2 = — 1. To prove this argument, let us perform the integration over |p| and |p'| in the 
small q limit, yielding 



/ 



d 2 pd 2 p' 4k 2 F f q 2 cos 6> pq cos fl p , q f pp , 

Cp/pp'Cp' = ——j / C)!t/p q dt/ p / q — — , (83) 



(2tt) 4 ^P p - p (2tt) 4 7 pq pq (n-v F q cos 9 p(l ){Sl - v F q cos ^ q ) 
where 6> pq (6> p / q ) is the angle between p and q (p' and q). In order to obtain the leading 
singularity near O = vpq, the above expression can be further simplified: 



f d z p d l p' 
J (27T) 4 C ^™' C *>' 



44 f M M /pp. 



de nn de 



(2tt)% 2 J ^ pq ^ p q [ (a_ _ A , 1^2 1 [fjL_iU 1^2 
For /pp/ oc ts — -75 — ix with A < 1, the above integral can be estimated through a scaling 

"pq W p'ql 



argument. We find 

d 2 p d 2 p' 



,.,_,( c p/pp' c p' 2+a 5 (85) 
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which is less divergent than (1 + a 2 ) 3 / 2 term that leads to 



Q 

v F q 



- 1 



-3/2 



divergence. 



Thus there is no cancellation between the first and the second terms in Eq.(82). If Svp 
diverges at small frequencies, we can conclude that 5il o will diverge in the limit v — > 
with u/vfQ fixed, which contradicts to our two-loop result from Eq.(71) that shows no 
such divergent term. Similar results also hold for the transverse current-current response 
function. 



The argument above assumes a power law behavior for / pp / oc 



Iflpq-VJ*' 



As A — > 1, 



another possibility needs to be considered, namely / pp / oc S(p — p'). This satisfies the 
condition that f 0s is finite. From Eq.(84) it is clear that this will lead to a term of 
order (1 + a 2 ) -3 / 2 which may cancel the first term in Eq.(82). However, in this case, 
we shall argue that, at least at zero temperature, / pp / = ( S(p — p') is equivalent to a 
shift in the Fermi velocity by vf — > v F + C^f/(27t) 2 . At zero temperature the excitation 
can be described by a distortion of the Fermi surface in the direction p by an amount 
Sup = f d\p\ Sn p . The original Landau's expression of the free energy density takes the 
form: 



SF = 



f d 2 p 1 f 

J J^)-2 M\P\ ~ k F )Sn p + - J 



(27T) 2 

Uf dp 1 

2 



d 2 p d 2 p' 

J PP 'dn p dn p/ 



2 v F (6vp) +- / -, , 4 fpp'dvpdvp 



(86) 



(2tt) 2 2 ry vj 2 J (2nY 
It is then clear that / pp / = C 5(p — p') is equivalent to vp — * vp + C^f/(27t) 2 . The same 
result can be also obtained by performing an integral over |p| in Eq.(79), which leads to 



(Q — vpq cos 9) Sup — q cos 9 



t/(q,0) + J fpp, Sup, 







(87) 



in the small q limit. Thus we see that, at zero temperature, all response functions to an 
external perturbation can be described by a Landau theory with a non-divergent effective 
mass in the small q limit. However, it is also possible that the same response function can 
be described by a Landau-Fermi-liquid theory of which both effective mass and / pp ' have 
divergent perturbative corrections. 
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An examination of Eq.(70) shows that after analytic continuation, the factor (1 + 
a 2 ) -5 / 2 diverges at O = vpq, even though its coefficient vanishes for fl — > 0. In the 
following we attempt an interpretation of the result. We can write our perturbative result 
Eq.(70) as, near fl = v F q, 

Tmnfnfll TmTT o (n I n dlmlQq^) i ,, a 2 Imn° (q,O) 

1m ll o(q, ") = lm ll 00 (q, S2) + a — h 70 , (88) 

where II[j is given by Eq.(77) with m* — > m, and 

a 2 

Q ' = ^' 



1 • // 1 11 f~<> j 2 



To = I I Q 

8n 2 (5 + r]) cog ^_2*_j k F ^ v F q \ \ ' 



4 



(89) 



where a 2 is a constant. The existence of d Im ITq (q, Q) /dQ term in Eq.(88) signifies 
that there is a finite non-singular (see a in Eq.(89)) shift in which also arises in the 
usual Fermi liquid theory. To interpret the second derivative term, we note that Eq.(88) 
is consistent with (apart from the term proportional to ao) 

Im n 00 (q, O) = 1 [ Im n° (q, Q + V) + Im Il^q, Q - T) ] (90) 

if T = v / 2to. We recall that Im n[] (q, ft) has a discontinuity at fl = VFq, corresponding 
to the edge of the particle- hole continuum. Eq.(90) has the natural interpretation of a 
smearing of the discontinuity at a shifted (due to a shift in vp) edge of the particle-hole 
continuum by the amount Y. Setting v F q oc O, we find that 

r oc ft 1+ ^ . (91) 

Note that for 1] < 3, V < fl so that the above picture is a self-consistent one. We also note 
that T is proportional to the square root of the coupling constant or 1/N, and is therefore 
non-analytic. We are not certain if any further physical meaning can be ascribed to the 
energy scale V. 
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VII. CONCLUSION 

In this paper we studied properties of gauge-invariant correlation functions in a two- 
dimensional fermion system coupled to a gauge field. We find the physical picture emerged 
from those gauge-invariant correlation functions to be very different from those obtained 
from gauge-dependent one-particle Green's function. The corrections to the Fermi-liquid 
two-particle correlation functions are found to be non-divergent and sub-leading to the 
Fermi-liquid contributions up to two-loop order, and there is no need to go beyond the 
perturbation theory at this order. 

However, it is still possible that singular corrections to the gauge-invariant two-particle 
correlation functions may appear in some special cases, such as q = 2kp- Also, since we 
do not have quasi-particles to serve as the underpinning of the Fermi-liquid-like behavior 
for n o and Iln, it is possible that singularity shows up in some other response func- 
tions. Nevertheless, the perturbative result should serve as a test for any theory such as 
renormalization group analysis [26] which attempts to go beyond perturbation theory. 

Finally we would like to comment on the implication of our results to the HTSC. 
Even though our results suggest that the two-particle Green's functions of fermions are 
Fermi-liquid-like for small q and O, it does not mean that the gauge field formulation of 
the t — J model (in relation to the normal state properties of HTSC) leads to the Fermi- 
liquid interpretation of the normal state of HTSC. In the problem of the t — J model, 
there are bosons as well as fermions which are interacting with a gauge field [12]. In fact, 
the presence of fermions and bosons in this problem came from the no-doulble-occupancy 
constraint on the electrons. It has been also regarded as a way of describing the spin- 
charge seperation induced by the strong correlation effects. In the paper of Nagaosa and 
Lee [12], they clearly demonstrated that the anomalous transport properties are due to 
the bosons. That is, the presence of the bosons plays an important role in the non-Fermi- 
liquid behaviors of the normal state of HTSC. However, in this paper we considered only 
the fermions interacting with a gauge field. 
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Figure captions 



Fig.l The one-loop diagrams for IIq (a) and for (b). The solid line is the bare electron 
propagator and the wavy line represents the gauge field propagator. These are the 
leading order diagrams of n o and IIn in the 1/N expansion. 

Fig. 2 The diagram that corresponds to the one- loop correction to the fermion self energy. 
The solid line is the bare electron propagator and the wavy line represents the gauge 
field propagator. 

Fig. 3 The diagrams that correspond to the (l/AT)°th order contributions to n n in the 1/N 
expansion. 

Fig. 4 The diagram that corresponds to the lowest order vertex correction r (k, q, iu, iv) or 
ri(k, q, iu, if). 

Fig. 5 (a) The non-vanishing diagram generated by ifj^a^A^ifj vertex, (b) A typical vanishing 
diagram generated by ifj'a^A^tf; vertex. 
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